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Abstract 

We explicitly evaluate the free energy of the random cluster model 
at its critical point for < q < 4 using an exact result due to Baxter, 
Temperley and Ashley. It is found that the resulting expression as- 
sumes a form which depends on whether 7r/2cos _1 (y / g/2) is a rational 
number, and if it is a rational number whether the denominator is an 
odd integer. Our consideration lealds to new summation identities 
and, for q = 2, a closed-form evaluation of the integral 
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1 Introduction 



The consideration of exact solutions of lattice models in statistical mechanics 
is a fertile field in both physics and mathematics. In physics it is well-known 
that (see, for example, pQ) exact solutions of lattice models lead to the finding 
of different classes of critical phenomena. In mathematics exact solutions 
in statistical mechanics often lead to the discovery of new mathematical 
identities. A prominent example is the solution of the hard hexagon model 
by Baxter [2] leads to the discovery of new Roger s-Ramanuj an summation 
identities in combinatorics [3]. As another example consider the integral 

I(A,B,C) = — I 2 * d9 d6\n\A + B + C 
An 2 Jo Jo 1 

-Acos6 - Bcoscf) - Ccos(9 + 0) (1) 

which arises often in lattice statistical studies. For v4 = I? = 2,C = 0itis 
well-known jH E] that the integral gives the per-site entropy of dimers on a 
square lattice and the integral can be directly evaluated |U Ej to yield 

s 4/111 \ 

For A = B = C = 2 the integral gives the per-site entropy of spanning trees 
on the triangular lattice, and the exact solution of the latter problem [HI Ej 
yields the closed form evaluation 

7(2,2,2) = ?Ii a (-L)+i In 3, (3) 

where Ti 2 is the inverse tangent integral function [S] defined by 

m , s f a tan" 1 t , 

Ti 2 (a) = / dt 

Jo t 

a 3 a 5 a 7 

Closed-form evaluation of I (A, B, C) for general A,B,C is not known. 

In this paper we report a closed-form evaluation of I (A, B, C), by using 
a result due to Baxter, Temperley and Ashley jH] on the random cluster 
model jUIl E] • We also show that the explicit expression of its critical free 
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energy depends crucially on the value of q. This implies that the free energy 
of the random cluster model, if solved, would also share this property. For 
convenience of references, we first state the integration result: 

For arbitrary A,B,C>0 and S = l/y/AB + BC + CA, we have 



I{A,B,C) = -ln(2S) + - 

71 



Ti 2 (AS) + Ti 2 (BS)+Ti 2 (CS) 



(5) 



2 The random cluster model 



Consider a g-state Potts model on the triangular lattice of N s sites with 
anisotropic (reduced) interactions K 1 ,K 2 , K 3 along the three principal lattice 
directions. Its partition function is 

Z^(K U K 2 , K 3 ) = J] HexpiKJ^Vi, a 3 )} (6) 

cr, = l E 

where the summation is over all q Ns spin states <7j, i = l,2,...,q, and the 
product is over nearest-neighbor set E connecting sites i and j, a = 1,2,3, 
and £>Kr is the Kronecker delta function. We are interested in the ferromag- 
netic regime K a > 0. 

The partition function can be rewritten jH] in a graph expansion as 

Z™(K u K 2 ,Kz) = Y / q <G \e K i - l)*(e*" - \)\e K * - 1)* (7) 

G 

where the summation is over all subgraphs G of the triangular lattice, c(G) is 
the number of connected clusters in G including isolated points, and £ a is the 
number of lines in G in the direction a = 1, 2, 3. In this form the partition 
function describes that of a random cluster model fDJ E] for which q can 
be continuous. One is interested in the per-site "free energy" of the random 
cluster model 

f RC = lim N; 1 In Z^(K U K 2 , K 3 ). (8) 
The triangular random cluster model is known to be critical at 

y/q Xix 2 x 3 + X\X 2 + x 2 x 3 + x 3 Xi = 1 (9) 

where 



Xr: 



-,K a 



i)/V?>o. (10) 
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For < q < 4, define 0(g) and v a (q), a = 1, 2, 3, by 

cos0(g) = y/q/2 , < < 7r/2 

x a = sin(0 — t> Q )/ smv a , < f Q < 0. 

For integral values of q, for example, we have 

q = 1, < w a < 0(1) 



9 
9 



2. 
3. 



< v a < 0(2) 



tt/3 
tt/4 



< v a < 0(3) = tt/6. 

In terms of these variables, we have 

1 



1 + 



g(4 — g) cot v a — q , 0<g<4 



(12) 



and the critical point © can be written as 

^1 + ^2 + ^3 = 20(g). (13) 

3 The critical free energy 

Baxter, Temperley and Ashley [Q^ considered the random cluster model (JTj) 
on the triangular lattice and obtained, among other things, its free energy at 
the critical point (JOJ). For < q < 4 they obtained 



J critical 



In q + ^(0, vi) + ip((f), v 2 ) + ip{(f>, v 3 ) 



where 



sinh(7r 



)x sinh 2(0 — v)x 



dx, < q < 4. 



(14) 



(15) 



xsinh7rx cosh (px 
Here, we investigate implications of ()14j) and (|15|). 

The integral (fTB^) can be evaluated by contour integration. Let C be the 
contour consisting of the real axis and the half circle at infinity encircling the 
upper half z plane. In the variable range < < 7r/2,0 < v < 0, we can 
neglect the contribution from the half-infinite circle and rewrite (|TH|l as 

)z sinh 2(0 — v)z 



if, n 1 f sinh(7r- 



zsmhitz cosh 02 



-dz. 



(16) 
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where the integration is alone the contour C. 

The contour integral can be carried out by using the residue theory. The 
integrand in (JTSJl has poles at 



sinh(7rzi) = 
cosh ((p z 2 ) = 



or Z\ = ni, 

71 

or z 2 = — (2m + lk, 
2q> 



n = 1,2, 



m = 0,1,2, 



The evaluation of the integral thus depends on whether z\ and z 2 overlap. 
Overlapping of z\ and z 2 occurs if 7r/20 is a rational number with the de- 
nominator an odd integer (M and N have no common integer factors): 



7T 



M 
iV' 



M = 1,2,3, 



N = 1,3,5, 



N < M. 



(17) 



Then there are two cases to consider. 



Case 1. There is no overlap between z\ and z 2 , namely, 7r/20 is either 
irrational or is of the form of (|T7|) but with iV even. This includes the q = 1 
Potts model (M = 3, N = 2). Then both z±,z 2 are simple poles, and the 
residues can be computed straightforwardly. One obtains 



ip((f),v) 



— tan(n0) sin 2n((f> — v) 



n=l 



n 



+ E 

m=0 



2m + 1 



• cot 



m + 



lx^ 



sm 



(2m + 1)— , (18) 



where the two terms come from the evaluation of residues at Z\ and z 2 , 
respectively. 

For later use the expression (|18|) can be further simplified by introducing 
in the first term the identity 

tanxsin2(a; — y) = cos 2y — cos2(x — y) + tanxsin2?/ 

with x = n<p, y = nv. This leads to 



ip((f),v) 



In 



sm 



v 



smv 

2 



OO J 

+ ^ — tan(n0) sin 2nu 



n=l 



r? 



+ E 

n=0 



2n + 1 



cot 



1X7T 2 



sm 



(2n + 1 



VTTi 



(19) 



after making use of the summation identity (1.441.2 of [TB*] ' 



E 

n=l 



cosnx 



111 



n 



2 sin (|) 



< x < 2ir. 



Case 2. There is overlapping between z\ and z 2 . This occurs when (II 7j) 
holds with iV = odd so that when 



2m + 1 = N, 3N, 5N ■ 



(20) 



we have z 2 = Mi, 3Mi, 5Mi, ■ ■ ■ which coincide with z\ where they form 
double poles. Note that when N — 1 every point in z 2 is a double pole. This 
includes the q = 2, 3 (M = 2, 3, N — 1) Potts models. The integral ([IB} now 
consists of three terms, ip 2 (NTr/2M,v) = Rsi + Rd + -R52, where Rsi is the 
sum of residues from simple poles in Zi, Rd the residues from double poles, 
and Rs2 the residues from simple poles in z 2 , if any. 

To compute Rsi, we note that the forming of double poles excludes points 
Mi, 3Mi, 5Mi, ■ ■ ■ which divide the remaining z\ = ni into sections n = 
{1, M - 1}, {M + 1, 3M - 1}, {3M + 1, 5M - 1}, ■ • ■. Then we can write 

Rsi = -Rsia + -Rsib (21) 

where Rsia is the sum of the first M — 1 residues and -Rsib is the sum of the 
rest. From we have already computed in the first term in (|18|). we obtain 

"sin[2/c(0 — v)]' 



R 



51a 



M-l 

tan(A;0) 

k=l 

M-l 00 



k 



R 



Sib 



V Y, tan[(2nM + k)<j>] 

-(Af-l) n=l 



^~2M 
sin[2(2nM + k)(<j> - v)] 



2nM + k 



Expanding the sine function, we can rewrite Rsi& as 



M-l 

#sia = E tan(fc^) 
fe=i 



/•n/2k r2v 

sm(2k(fi) / sin fcxdx — cos(2&;</>) / cos /cxiia; 
J2d Jo 



To evaluate -Rsib we use 4M0 = 2iW and proceed similarly as in 
obtain 



(22) 
and 



M—l 00 

Rsib = tan(A;0) ^ 

fc=-(Af-l) n=l 



sin(2A;0) 



7r/2|fc| 



2 c 



/■2w 

-cos(2A;0) / cos[(2nM + fc)a;]dx 
Jo 



sin[(2nM + fc)x]daj 
, Nit 



2M 



where one verifies that the integrated contributions at the upper limit n/2\k\ 
cancel out exactly after combining terms from k and —k. Here the k = 
term vanishes due to the tan(fc0) factor. 

Expanding out the sine and cosine functions as before, terms odd in k in 
the summand are canceled and terms even in k give the same contribution 
for k > and k < 0. Thus we have 



A/-1 oo 

Rsi b = 2 tan(fc0) £ 

k=l n=l 



sin(2fc0) 



7r/2fc 



2t' 



sin(2nMx) cos(kx)dx 
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Write 



+ cos(2k(f)) / sin(2nMx) sin(kx)dx 



oo Af 

sin(2nMx) = lim ^ sin(2nMa;) 

n=l °° n=l 



(23) 



lim 

Af— >oo 



cos(Mx) - cos[(27V+ l)Mx] 
2 sin Mx 



(24) 



where we have used the summation identity 1.342.1 of jTSj hi writing down the 
last step. Substituting into (J23j) . terms with the factor cos[(27V+ l)Mx] 
inside an integral vanish in the J\f — > oo limit. Thus, one obtains 



M-l 

-Rsib = XI tan ( fc 0) 
k=i 



sin(2A;0) 



vr/2fc 



2 u 



cos(/ca;) cot(Mx)dx 
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Finally, combining 



M-l 

R si = tan(fc^) 
fc=i 



+ cos(2A;0) / sin(fcx) cot(Mx)dx 
Jo 

and (jUJ), we obtain 

^/ 2fc cos(M - k)x 



(25) 



sin(2/c0) 
— cos(2A;0) 



2d sin(Mx) 
2u sin(M - k)x 
o sin(Mx) 



-dx 



dx 



Ntt 
2M' 



(26) 



Now we compute Rd- Residues from double poles at z 2 = (2m + l)Mi 



are 



00 d 
Rd = 7ii Y Ihri — 

* — ' 2^22 (17 

m=0 UZ 



yZ - Z 2 



( sinh(7r — (j))z sinh 2( 



VIZ 



zsinhvrz cosh chz 



(27) 



where z 2 = (2m + l)Mi, = Nn/2M, N = odd. 
Define number u and integer p by 

Mv=pn/2 + u, where < u < 7r/2, p = 0, 1, 2, • • • , p < N. (28) 

Here for technical reasons we exclude points u = 0. At these points each of 
Rgi, Rd and Rs2 diverge. But in this cse -0(0, v) can be computed by taking 
an appropriate large M, N limit of (fT^j). (See discussions in Section 5 below.) 
Then using the identities 

Km 4^ = -("1) M lim = ^(-l)™+<"-M» 

sinh 7rz 7r 2 ^ z 2 cosh 0z AT 7T2 

sinh(7r-0)z 2 = z (_i)M-m+(Af+i)/2 ; smn 2 (0 - t;),2 2 = z (- 1 ) p sin 2 (2m + 1 )u 
cosh(vr - 0)z 2 = 0, cosh 2(0- v)z 2 = (-l) p+1 cos2(2m + l)u 

it is straightforward to obtain from (}2*Tj) 

_ 2(-l) p ~ sin[2(2m + l)u] 4(-l) p (0 - v) ~ cos[2(2m + ]>] 
D "MA%^ (2m +1)2 + A% ^ 2^Tl ' 

(29) 

The first summation in (|2*9*|) can be carried out by using the summation 
identity 3.1 established below, and the second summation carried out using 
1.442.2 of H3|, namely, 

^ cos 2(2*;+ l)x 1, 7T 

Y )■ '— = - In cots, 0<a;<-. 30 

to 2k + l 2 2 1 J 

This leads to the simple expression 

2(-l) p , v (-l)P(N-p), , . 



Identity 3.1. 



~ sin2(2A; + lk m , vr . . 

2 y — t— : -Tj — = Ti 2 (tanx) + xlncoti, < a; < — . (32) 

fc =o ( 2fc + 1) 2 



Proof. We have 



~ sin 2(2A; + l)x ~ cos(2fc + l)y 



fc=0 (2A;+1)2 7o fc f 2fc + 1 

1 ,3* , y 



In f cot — )cfy 
2 7o v 2^ y 

/■tan x ] n + 

/ dt , < x < -, (33) 

io l + t 2 ' 2' v ; 
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where we have made use of (|3Uj) and the change of variable t = tan(y/2). 
This leads to after using the identity (4.29) of [B\, namely, 

rt&nx ] n f Tj- 

Ti 2 (tanx) = xln(tana;) — J - — ^ dt, x < — . (34) 

Q.E.D. 

To compute R$2, the sum of residues of simple poles in z 2 , we need to 
exclude from z 2 — [(2m+l)M/N]i the double poles Mi, 3Mi, 5Mi, • - -. As a 
result, the remaining simple poles in z 2 are divided into sections m = {0, (iV— 
3)/2}, {(N+ l)/2, (3iV - 3)/2}, {(3iV + l)/2, (5iV - 3)/2}, ■ ■ -. The situation 
is similar to that of Rsi, and Rg2 can be similarly computed. Omitting the 
details we arrive at the expression 

2M (A ^ )/2 ^ (2kMix\ f*> sm(2kMx/N) 

Rs2 - 



N k 



~" /2kMTT\ r M sin [2k Mx UN , , . 

E cot (— aT" / • 7T7 ^ V iV=l,3,5,---.35 

\ N J Jo sm(Mi) 



Note that i?s 2 = when AT = 1, as every point in z 2 is a double pole. 
Finally, combining the above results, we obtain the expression 

K^ ,V ) = R si + Rd + R S 2, M = l,2,3,--- AT=1,3,5,--- 



where Rsi is given in (J2EJ)> Rd by (}3*T]) . and _R 52 by (|3*5]). 



(36) 



4 Special cases 

In this section we consider the special cases of q = 1, 2, 3 Potts models. 

For q — 1, we have = 7r/3, M = 3, A" = 2, so we use (fTTIj) . For 
= 7r/3, we have cot [(n + l/2)7r 2 /0] = and the third term in (j!9|) vanishes 
identically. It is then a simple matter to use the summation identity 4.1 we 
establish below to obtain 

»(f,«)=h ^' + 1 . (37) 
It follows that the free energy is 

f RC = K 1 + K 2 + K 3 , q = l 
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after making use of (|12j) and (|14|) . This result agrees with the physical 
expectation that, since there is only one spin state in the q = 1 Potts model, 
the partition function is Z = exp[N s (Ki + K 2 + K 3 )]. 

Identity J^.l. 

cot u + 1 , , s 

> — tanm7r/3) sin 2m; = In — = , < v < 7T/6. (38) 

^ n v 3 cot u — 1 

Proof. Since tan(rar/3) = for n = 3 x integers, we have 



n=l 



— tan(n7r/3) sin 

n=0 

00 

= v^E 



sin 2(3n + l)v sin 2(3n + 2)u 



3n + 1 



3ri + 2 



n=0 



2(! 



cos(3n + l)x — cos(3n + 2)x 



dx 



v3 lim V" / cos(3fc + l)x — cos(3/c + 2)x 



k=0 



dx 



The summation can be carried out by using the formula (1.341.3 of |13j ) 



n-l 



v-» / , , \ ( n-l \ . ny y 
} t cos(x + ky) = cos x H — y I sin — esc - . 

k=o \ 2 / 2 2 



This yields 

00 1 



E — tan(rar/3) sin2nt> = v3 lim / 
~[ n JV-+00 7o 

= Vs 



[ 1 — cos 3J\fx) sin(x /2) 



2v 



2 cos a; + 



sin(3x/2) 
y/3 lim / 

1 Af^oo Jq 



dx 

cos 3Afx 
2 cos x + 1 



dx. 



(39) 



The second term in (J59j) vanishes and the first term can be evaluated by 
using the formula (Cf. 2.551.3 of [T3j ) 



dx 



1 



a + 6 cos x a/6 2 — a 2 
This leads to (jMD- Q.E.D. 



In 



V& 2 — a 2 tan(f) + a + b 
\/b 2 — a 2 tan(|) — a — b 



a 2 <b 2 
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For q = 2 we have = 7r/4, M = 2, N = 1. By comparing our result 
with that of the Ising model, the critical free energy leads to the integration 
formula (jSJ) as we shall now show. 

First, using we obtain 

ip(^,v) = -]n[(cot«)(cat2u)] + -Ti 2 (tan2v), q = 2. (40) 

\4 / 2 7T 



Therefore from (|14jl the critical free energy is 
1 3 



1 a=l 



^ln[(cotu a )(cot 2v a )\ + -Ti 2 (tan2f a ) 

2 7T 



9 = 2. (41) 



On the other hand, the q = 2 Potts model is completely equivalent to an 
Ising model. Consider an Ising model on the same triangular lattice of N s 
site with anisotropic interactions K a /2, a = 1, 2, 3. We have the equivalence 



Zn7 9 = e n 

<7=±1 E 



D (K a /2)tTitTj 



e -N s (K 1 +K 2 +K 3 )/2 ^Potts 



g=2 



(42) 



after making use of the identity a^j = 25 K r(c r i, &j) — 1- It follows that their 
critical free energies are related by 



/ 



Ising 



f 



Potts 



f Potts\ i n [( cot Ul )( (5 ot u 2 )(cot v 3 )] 

I q — 2 2i 

i 3 r l l 

-hi2 + ^ -ln(cot2t; a ) + -Ti 2 (tan2< 



after substituting with (j4~Tj) and using (|T2*j) for g = 2, namely, 
e x " = cot u Q , or sinh K a = cot 2iv 



(43) 



(44) 



In fact, the Ising free energy is known at all temperatures ^3] to be 



/ 



Ising 



In 2 + — / / In 



cosh ^ cosh K 2 cosh + sinh K\ sinh sinh K 3 

(45) 



87?' JO JO 

— sinh Ki cos # — sinh K 2 cos — sinh K 3 cos(8 + 
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It can be verified that the critical condition JSJ) is equivalent to 

cosh K\ cosh K 2 cosh K$ + sinh K\ sinh K 2 sinh K% 

= sinh K\ + sinh K 2 + sinh i^, (46) 

so we can rewrite (|45|) as 

fising = i n2 + -I(a,b,c) (47) 

with a = cot 2vx, b = cot 2t> 2 , c = cot 2w 3 and /(a, 6, c) defined in (JTJ . It can 
also be verified that the critical condition v 1 + v 2 + V3 = 7r/2 is the same as 



ab + be + ca = 1 . 
Equating (J4l?|) with (J4TJ) . we obtain 



(48) 



I{a,b,c) = -In2 + ln(a6c) + - Ti 2 (a _1 ) + Ti^Zr 1 ) + Ti 2 (c 



7T 



(49) 



For the general integral I(A, B, C) we can always define variables a = 
AS,b = BS,c = CS with S = 1/y/AB + BC + CA to satisfy flUD- Thus, 
one has 



I(A,B,C) 



In S + I (a, b, c) 



2 r 



= - \n(2S) + m(afec) + - Ti^a" 1 ) + Ti^ir 1 ) + Ti 2 (c 

7T L 

= -ln(2^) + -[Ti 2 (a)+Ti 2 (fo)+Ti 2 (c)l, (50) 

where use has been made of the identity (Cf. (2.6) of |H]) 

Ti 2 (y- l )=Ti 2 (y)-^\ny, y > 0. 

The last line in (|50)l establishes the integration formula JHJ). It is readily 
checked that (JHOj) yields the previous known values of 7(2, 2, 0) and 1(2, 2, 2). 

For q = 3 we have = 7r/6, M = 3, N — 1, so we again use (JHEJ- It is 
straightforward to obtain 

v^cot 2v + 1 



' '( 



1, 2cos2w + v / 3 3, v^+l 
- m p= H — In 



In 



6 2cos2w- v / 3 2 y/3 - 1 V / 3cot2w-l 

2 \/3 12 7r 

+ — In 7= H — lncot3iH Ti 2 (tan3t>), v < — . (51) 

12 2 + V3 3 37r 6 
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5 Summary and acknowledgments 



We have evaluated explicitly the free energy of the random cluster model 
at its critical point for < q < 4 using an exact result of Baxter, Temper- 
ley and Ashley. It is found that the critical free energy is given by (fTK|) if 
7r/20(g) = 7r/2 cos _1 (y / g/2) is irrational or a rational number with an even 
denominator, and by (j3T?j) if 7r/20 is a rational number with an odd integer in 
the denominator. Special cases of our results lead to new summation identi- 
ties (J3~2|) and (pl%j) . and a new integration formula © for I (A, B,C), which 
do not appear to have previously appeared in print. 

It is instructive to see how the critical free energy passes from (|3l)|) to 
(|TH|) as q varies and 7r/20 changes from rational to irrational. Indeed, any 
irrational 7r/20 can be reached by taking an appropriate M, N — > oo limit of 
7r/20 = N/M. In this limit we have Rp = by 1)31)1 . It can be verified that 
Rsi and R$2 in and (J3~5j) can be rewritten in equivalent forms 



Af-l 



R 



SI 



— tan(n0) sin[2n(0 — v )] 



n=l 



n 



Af-l 



+ ^ tan( 



n=l 



7r/2n 

sin(2n0) / cot(Ma;) cos(na;)(ia; 



2r 



2 c 



+ cos(2n0) / cot(Mx) sm(nx)dx 
Jo 



R 



52 



(N-3)/2 

2 51 cot 

m=0 



m + 



1\ 7TX 



sin[(2m + l)vn/<p] 



2m + 1 



AT 7 



cot (Ma;) sin 



m + 



l\nx 



dx 



(52) 



In the large M, N limit, the integrals in ()52)l vanish. Then R$i becomes the 
first term in 1)18)1 . i?^2 becomes the second term in ()18jl and one recovers l)T8)) 
from ()3f)jl . Similarly (p?fi)) can be reached by taking an appropriate limit of 
fill?)) or (JTHJ. Particularly, in the case of the points u = excluded in (J2BJ)> 
i/j((f),v) can be compute by taking an appropriate large M,N limit of ()19p. 

We would like to thank Professor Shin-Nan Yang for the hospitality at 
the Center for Theoretical Physics, Taipei, where this work is initiated. LCC 
is supported by a travel grant of the Institute of Mathematics, Academia 
Sinica, Taipei. FYW wishes to thank W. T. Lu for a useful conversation. 
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